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1. Introduction 

The properties of interfaces in three-dimensional statistical systems have been a long-standing 
subject of research. In particular the interest of people working in the subject has been attracted 
by the so called "fluid" interfaces whose dynamics is dominated by massless excitations (for a 
review see for instance [1]). For this class of interfaces, thanks to the presence of long range 
massless modes, microscopic details such as the lattice structure of the spin model or the chemical 
composition of the components of the binary mixture become irrelevant and the physics can be 
rather accurately described by field theoretic methods. 

A simple realization of these fluid interfaces is represented by 3d spin models. These models 
in the low temperature phase in which their global symmetry is broken admit different vacua which, 
for a suitable choice of the boundary conditions, can occupy macroscopic regions and are separated 
by domain walls which behave as interfaces. For temperatures between the roughening and the 
critical one, interfaces are dominated by long wavelength fluctuations (i.e. they exactly behave as 
fluid interfaces). 

An effective model widely used to describe a rough interface is the capillary wave model 
(CWM) [2, 3]. Actually this model (which was proposed well before the Nambu-Goto papers [4, 5]) 
exactly coincides with the Nambu-Goto one, since it assumes an effective Hamiltonian proportional 
to the variation of the surface's area with respect to the classical solution. 

Among the various realizations of fluid interfaces, a prominent role has been played in these 
last years by the 3d Ising model, for several reasons. The universality class of the Ising model 
includes many physical systems, ranging from binary mixtures to amphiphilic membranes. It also 
shares the universality class with the 0^ theory; this allows a QFT approach to the description of 
the interface physics. Last but not least, the Ising model, due to its intrinsic simplicity, allows fast 
and high statistics Monte Carlo simulations, so that very precise comparisons can be made between 
theoretical predictions and numerical results. 

Another important reason of interest of interfaces in the 3d Ising model is that, thanks to the 
duality relation with the 3d gauge Ising model, they are deeply related with Wilson loop observ- 
ables. This is the ultimate reason for the validity of the Nambu-Goto effective string model for 
interfaces. Indeed this mapping allows to create a dictionary between Lattice Gauge Theory (LGT) 
observables and interface ones and makes the results obtained in the context of interface physics 
immediately relevant also for LGT studies. In this respect, it is important to stress that by working 
directly in the 3d Ising model one has a larger freedom in the choice of the boundary conditions. In 
this work we shall exploit this freedom and address interfaces with periodic boundary conditions. 
This choice is particularly useful if one wants to address subtle finite size connections (as one has to 
do in order to test the Nambu-Goto effective string picture) since the periodic b.c. allow to elimi- 
nate any source of boundary corrections in the results, thus leading to a very clean and artefact-free 
comparison between theory and simulations. 

The aim of this contribution is to summarize a set of results recently obtained in this context 
by our group [6, 7, 8, 9]. In particular we shall discuss the following items: 

1] Explicit evaluation of the Nambu-Goto partition function (neglecting the anomaly) for inter- 
faces [6]. 



2 



Effective string theory description of the interface free energy 



Michele Caselle 



2] A collection of numerical results from a set of high precision Montecarlo simulations [7, 8] 
3] The 2d limit of 3d interfaces [9]. 

2. The partition function of interfaces from the Nambu-Goto effective string 

The Nambu-Goto bosonic string is defined by an action proportional, via the string tension a, 
to the induced area of a surface embedded in a J-dimensional target space: 

f ^ r. dX' dXj 

S = ajd-^y^g, g„p = ——Gij, (2.1) 

where the proper coordinates parametrize the world sheet while X'{^) (M = 0,.. .,d) describes 
the target space position of a point specified by ^ . We shall assume in the following that the target 
space metric Gjj is the flat one and neglect it. 

The standard approach to evaluate the partition function is to fix the re-parametrizations of the 
action by fixing the so-called "physical" gauge in which the proper coordinates are identified with 
two of the target space coordinates, say X^ and X'. The quantum version of the NG theory can 
then be defined through the functional integration over the J - 2 transverse d.o.f. X{X^\X^) of the 
gauge-fixed action. The partition function for a surface L with prescribed boundary conditions is 
given by 



Zi: = J^^ DX' exp | -(^^ dX^dX^ (^1 + (^o^)^ + {diXf + {doX A dyX) 

= [ DX' expi-a / dX^dX^ 
J{dz) [ Jz 



1 + ^ {d^Xf + ^ {diXf + interactions 



(2.2) 



Expanding the squai^e root as in the second line above, the classical area law exp(— aiz/) 
(where ^ is the area of the minimal surface £) is singled out. It multiplies the quantum fluctua- 
tions of the fields X, which have a series of higher order (derivative) interactions. The functional 
integration can be performed perturbatively (since we are studying an effective model we can ne- 
glect the fact that the theory is actually non-renormalizable), the loop expansion parameter being 
1 /(a^), and it depends on the boundary conditions imposed on the fields X, i.e., on the topology 
of the boundary dlL and hence of Z. The cases in which £ is a disk, a cylinder or a torus are the ones 
relevant for an effective string description of, respectively, Wilson loops, Polyakov loop correlators 
and interfaces in a compact target space. The computation was carried out up to two loops in [10]. 

The main problem of this approach is that the physical gauge is anomalous in J / 26 dimen- 
sions. However it can be shown [11] that the anomaly decreases as R^, where R is the typical 
scale of the boundary of £ and thus it could be possible that it does not affect the first two terms in 
the loop expansion. 

In [6] we proposed an alternative treatment of the NG model which takes advantage of the first 
order formulation, in which the action is simply 

S = o j d^"^ j^"" d^^h^^daX'dpX' , (2.3) 



3 



Effective string theory description of the interface free energy 



Michele Caselle 



where /j„jg is an independent world-sheet metric, G [0,27r] parametrizes the spatial extension 
of the string and i^" its proper time evolution. Integrating out h, we retrieve the NG action eq. 
(2.1). For each topology of the world-sheet, classified simply by its genus g, we can instead use 
re-parametrization and Weyl invariance to put the metric in a reference form e'^ha^ (conformal 
gauge fixing). For instance, on the sphere, i.e. at genus g = 0, we can choose h„j} = rjap, while on 
the torus, at genus g = I, hap is constant, but still depends on a single complex parameter T, the 
modulus of the torus. The scale factor e*^ decouples at the classical level and the action takes then 
the form 



S = O J J^'' d^^ h^^daX'dpr + 5gh. , (2.4) 

where ^gh. in eq. (2.3) is the action for the ghost and anti-ghost fields (traditionally called c and 
b) that arise from the Jacobian to fix the conformal gauge; we do not really need here its explicit 
expression, the only important thing to know is that they correspond to a CFT of central charge 
Cgh. = —26. The fields Z'(t, a), with i = \, . . . ,d, describe the embedding of the string world-sheet 
in the target space and form the simple, well-known two-dimensional CFT of d free bosons. This 
allows to evaluate exactly the sum over all possible surfaces for any value of d. In the particular 
case of the toroidal geometry (interfaces) in which we are interested one finds (see [6] for the 
details): 



j{d) = 2 ( ' Vt Vas^u £ CkCk' ( - ) " Kd^ {a£/E) , (2.5) 



2nJ 



k.k'=0 



d-l 
2 



where .s^ = L\L2 denotes the minimal area of the interface, m = ^ is the ratio between the two 
lattice sizes of the plane in which the interface lies (L, denotes the lattice sizes in the /th direction), 
Vt = nf=3 Li is the transverse volume, Ck is the ^th coefficient in the expansion of the Dedekind 
function: 



[n{qf-'=Y,c,q'-'^ , (2.6) 

k=0 

Kd-i is the modified Bessel function of order and 

2 

denotes the spectrum of the NG string for this particular type of geometry. 

It is important to stress that also in this case for J / 26 we must expect deviations from this 
result. In fact the decoupling of the scale factor e'^ persists at the quantum level only if the anomaly 
parametrized by the total central charge c = d — 26 vanishes; Weyl invariance is otherwise broken 
and the mode (p has to be thought of as a field with a Liouville-type action. 

Notwithstanding this there are a few non-trivial results of this analysis which are not affected 
by the anomaly and deserve further attention: 

• Eq.(2.5) gives an explicit prediction for the functional dependence on a, £/, u and Vt of the 
prefactor in the interface free energy which (as we shall see in the next section) nicely agrees 
with the results of the MC simulations. 
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square interfaces Lq=96, L2=64 

P=0.223102 




Figure 1: Fig.la: Comparison of the 2-loop prediction, the full Nambu-Goto result with the data for square 
interfaces at j3 = 0.223102 and j3 ^ 0.226102. Fig.lb: F,'^' - (7L1L2 as a function of y/aLi for )3 = 
0.223102, L3 — 96 and L2 — 64. Note that in the case of the Monte Carlo results the statistical eiTor is 
smaller than the symbol (circle). The 1-loop, 2-loop and full NG predictions are given as solid black, red 
and blue lines, respectively. 



• The spectrum reported in eq.(2.7) is different from the one obtained by Ai vis in [12] for the 
Wilson loop geometry and agrees with an independent calculation mentioned in [13]. 

• Eq.(2.5) in the d = 3 case allows a smooth dimensional reduction hmit and suggest an ex- 
pression for the 2d interface which agrees with an exact 2d calculation. 

• Expanding eq.(2.5) in powers of 1/{G£/) we reobtain at the first two orders the same con- 
tributions obtained in [10] with the physical gauge. This is probably the most interesting 
result since it supports the idea that not only the first term in the expansion (the so called 
"Liischer term" [14]) but also the next to leading order could be anomaly free and universal. 
This conjecture is further supported by the recent observation [15, 16] that the same term 
also appears in the expansion of the Polchinski-Strominger action [17] which is essentially a 
non-local, but anomaly free, version of the N-G action. 

3. Monte Carlo simulations of interfaces in the three-dimensional Ising model 

In order to test the above results it is mandatory to have very precise estimates of the interface 
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free energy for different values of Li and L2. In [7] and [8] we obtained these estimates using two 
different approaches. In [7] we used the "boundary flip" algorithm [18] which is particulai^ly suited 
for interfaces of small or intermediate size while in [8] we obtained the interface free energy by 
a direct integration of the differences — Ep over j3 (where Ea denotes the internal energy for 
antiperiodic boundary conditions, and hence in presence of an odd number of interfaces, while Ep 
denotes the internal energy with periodic b.c). This second method is much more cumbersome 
from a computational point of view, but it is mandatory if one is interested in large values of the 
interface free energy. We refer to the original papers for all the details on the algorithms and on the 
simulation settings and only list here the main results. 

• As it was already observed in the case of Wilson loop [19] and of Polyakov loop con^ela- 
tors [20], we found a very good agreement for large values of R and a sizeable disagreement 
as R decreases (see fig. la). This deviation (due to the peculiar geometry of the interface 
observable) cannot be due to boundary corrections and is thus an intrinsic feature of the 
effective string description. In the case of square interfaces (i.e. L\ = L2 = R) the NG ef- 
fective string picture breaks down for values of R smaller than R ~ 2.5 / \/a (see fig. la). It 
is important to notice that for square interfaces we see essentially no difference between the 
prediction obtained using eq.(2.5) (black line in fig. la) and the two loop approximation (red 
line in fig. la). 

• In order to disentangle between the two predictions one must study asymmetric interfaces 
because higher orders in eq.(2.5) become more and more relevant as w = L2/L1 increases. 
This is evident in fig. lb where the two curves are well separated. Looking at fig. lb we 
see that the MC data show a much better agreement with the two loop approximation than 
with eq.(2.5). This deviation from the behaviour of eq.(2.5 may simply denote a breaking 
of the effective string picture in this asymmetric limit. However it is also possible (and the 
persistent agreement with the two loop approximation could support this interpretation) that 
the effective string picture still holds but the corrections due to the anomaly (or equivalently, 
in the covariant gauge framework, to the Liouville field) which affect the calculation of 
eq.(2.5) cannot be neglected. 

• Last, but not least, thanks to the improvement in the estimate of the interface string tension 
a we could obtain new predictions for a few amplitude ratios which are more precise than 
any other existing estimate in the literature: 




1.235(2) 



mo++ 



3.037(16) 



R-^ 0^1, = 0.1024(5) 
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4. Universal behaviour of interfaces in 2d and dimensional reduction of 
Nambu-Goto string 

It is interesting to study the m — > oo (i.e. Li << L2) limit of the eq.(2.5) for d = 3. This 
corresponds to perform a dimensional reduction from d = 3 to d = 2 and allows us to obtain an 
effective description for bidimensional interfaces. The result is [9]: 



Eq.(4. 1) is most probably an exact result for any 2d spin model and not just an effective approxi- 
mation. Indeed it agrees with an exact calculation in the 2d Ising model. Moreover it arises from 
a simple universal model for 2d interfaces, based on the first order treatment of an action given by 
the length of the interface. 

Acknowledgements: M.P. acknowledges support from the Alexander von Humboldt Founda- 
tion. 
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